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Abstract 

RGIsearch is a C++ program that searches for invariants of a user-dehned 
set of renormalization group equations. Based on the general shape of the 
+functions of quantum held theories, RGIsearch searches for several types 
of invariants that require different methods. Additionally, it supports the 
computation of invariants up to two-loop level. A manual for the program 
is given, including the settings and set-up of the program, as well as a test 
case. 

Keywords: Renormalization, Computer Algebra, Sparse Linear Systems 


* Corresponding author 

Email address: RobVerheyen@gmail.com (Rob Verheyen) 

BMAPP, Mailbox 79, Radboud University Nijmegen, P.O. Box 9010, 6500 GL 
Nijmegen, The Netherlands 


Preprint submitted to Computer Physics Communication 


July 15, 2015 





Program Summary Manuscript Title: RGIsearch: A C++ program for 
the determination of Renormalization Group Invariants 
Authors: Rob Verheyen 
Program Title: RGIsearch 
Journal Reference: 

Catalogue identifier: 

Licensing provisions: none 
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Operating system: Linux 
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Nature of problem: The determination of renormalization group invariants, which 
can be used as a probe for the high-energy behaviour of theories in particle physics. 

Solution method: Several types of invariants are considered based on the general 
shape of the renormalization group equations. The problem of computing these 
invariants is then reduced to creating and solving very large systems of linear 
equations or generalized eigenvalue problems. Using the specific sparsity structure 
of these systems, the systems can be solved. 

Restrietions: Since the algorithms that solve the linear systems and generalized 
eigenvalue problems are very efficient, the main restriction is the amount of 
available RAM, where the systems are stored. The program supports polynomial 
renormalization group equations, which are typical particle physics, up to two-loop 
order. 

Running time: Ranging from seconds to minutes. 
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1. Introduction 

Modern theories in particle physics, such as supersymmetry, often predict 
new physics at experimentally inaccessible energy scales. The high energy 
behavior of such a theory can be probed by considering the evolution 
of its parameters, which were measured at some experimental scale, to 
higher energies using the renormalization group (RG) equations. Low-energy 
parameters can be translated to high-energy ones with several methods, each 
with their own up- and downsides [ij. One of these methods uses RG 
invariants, combinations of the parameters of a theory that are invariant 
under the renormalization group flow While this method offers multiple 

advantages over others, hnding such invariants is highly nontrivial, especially 
at high loop orders or for theories with large parameter spaces. RGIsearch 
is a G-I--I- program that is able to compute these invariants for any set of 
RG equations. In this article, the methods that RGIsearch uses to hnd 
invariants are described and a manual is provided. Obviously, checking if 
the proposed invariants produced by the program are actually invariant is 
trivial. RGIsearch was recently used to compute several new invariants of 
various versions of the Minimal Supersymmetric Standard Mode j^. 

2. Methodology 

The RG equations are differential equations for the parameters of a theory 
as function of the energy scale /i. Consider the RG equation of a renormalized 
theory for a parameter x{fi). The corresponding /3-function is dehned as: 



( 1 ) 


where t = log^g (h/ho)- Since the /3-functions of quantum held theories 
are polynomials in the parameters of the theory, RGIsearch considers several 
types of polynomial RG invariants. We hrst discuss methods for these types 
for one-loop /3-functions, where in quantum held theories the coefhcients of 
these polynomials are rational numbers, before extending it to two loops. 

2.1. Monomial Invariants 

Monomial invariants assume the form: 


n 


M=n 



( 2 ) 


2 = 1 
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where Xi are the n parameters of the theory and a E 7/^. The powers a can 
be taken to be integers since the coefficients of the /^-functions are rational, 
and any power of a monomial invariant is still a monomial invariant. Taking 
the derivative and dividing out M, the condition for invariance is: 

= 0. (3) 

• 1 
1=1 

The left-hand side of eq. (jS]) is a polynomial which is required to vanish for 
all values of all Xi. This is only possible if all monomial terms are cancelled 
internally. Therefore, eq. ([2]) can be reduced to a relatively small linear 
system of equations on the powers a*, where the equations represent the 
requirement of the cancellation of all monomial terms. Every element in the 
nullspace of this linear system represents a monomial invariant. By Ending 
the basis vectors for the nullspace, the independent monomial invariants are 
determined. The other elements of the nullspace are then just products of 
powers of these invariants. 

2.2. Dimensionalities 

Before extending the above method to polynomial invariants, the concept 
of dimensionalities has to be introduced. A dimensionality D is represented 
by a vector D E 7'^ such that: 


dimoixi) = Di. 

Dimensionalities are additive for monomials: 


( 4 ) 


diniz) ( =a- D, 


( 5 ) 


V 2=1 


where • represents the standard Euclidean inner product. A system of 
/9-functions is said to have dimensionality D if: 


Vz e [1, ...,n] : dimB(/3(a:j)) - dim£,(a;j) = cd, (6) 

where the dimensionality of a polynomial such as the /9-functions means 
that all included monomial terms have the same dimensionality, and cd E 7 
is a constant. A system of /9-functions can have multiple dimensionalities. 
Note that these dimensionalities are similar to the well-known physical 
dimensionalities, such as mass dimension. In fact, mass dimension is typically 
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one of the dimensionalities of a system of /^-functions with = 0. 

RGIsearch hnds the dimensionalities of a system of /9-functions by reducing 
the problem to a different system of linear equations. Writing out the 13- 
functions in their monomial terms: 

m.i n 

where (7) 

j=l k=l 

where the rrii are the number of monomial terms in f3{xi). Eq. ([6]) can 
now be written as: 


Vi e [1, n], j e [1, rrii] : kj ■ d - Di - Cd = 0, (8) 

which is a linear system in d and cd and can thus easily be solved. 


2.3. Polynomial Invariants 

Similar to the /9-functions, polynomial invariants assume the form: 


m n 

P = ^ CiMi where > 

i=l j=l 

where C G Z'^ and a, G Z". Eq. ([6]) implies: 


( 9 ) 


dim^) 



dirni) (Mi) -1- Cd- 


( 10 ) 


This means that there cannot be any overlap between derivatives of 
monomial terms with different dimensionalities. Therefore, it suffices 
to consider polynomial invariants with monomial terms of the same 
dimensionalities. Thus, RGIsearch first computes the dimensionalities of the 
system. Next, if r dimensionalities are found, it runs through a user-defined 
range of dimensionalities d G Z*, computing the set: 


Mp(d) = }^W_x'l^\\/le[l,...,r]-. dim, ( j = dz, -p < Xi <p I , (11) 


. ^ = 1 


v2=l 


where p is a user-dehned setting controlling the size of A4p(d). The 
polynomial invariant then becomes: 
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( 12 ) 




P{d) = CiMi where Mi G A4p(d). 


i=l 


The requirement for invariance now is: 



i=l 


(13) 


which can be converted to a linear system of equations in C by demanding 


that eq. flT^ holds for all values of all parameters. While the systems 
encountered previously were relatively small, the size of this system heavily 
depends on the size of Mp, which should be taken as large as possible to 
reach more invariants. The associated system of equations can grow very 
large, but it can still be solved quickly due to its sparsity. 

The elements of the nullspace make up all possible invariants that can 
be constructed from the set Mp. By computing the basis vectors of the 
nullspace, RGlsearch hnds all independent invariants. The rest of the 
nullspace consists of linear combinations of these invariants. In addition, 
products of these invariants can be found at other dimensionalities. After 
running through the range of user-dehned dimensionalities, a simple hltering 
algorithm searches for these products of invariants and removes them from 
the output. 

S.4- Factorization 

RGlsearch supports a minor extension to the method of the previous 
section, where a single variable can be factorized from a polynomial invariant: 



Taking the derivative and dividing out the requirement for invariance 
is: 



( 15 ) 


Instead of a linear system, eq. (irsj) leads to a system of the form: 
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A + feQ = 0. (16) 

This can be interpreted as a nonsqnare, generalized eigenvalne problem of 
a similar size as regnlar polynomial searching. It is approached as a regular 
eigenvalue system by solving the eigenvalues b hrst, after which the system 
becomes linear. 

2.5. System Creation 

The methods described in subsections 2.3 and 2.4 require the creation of 
very large (typically 0(10^) x 0(10^) to (9(10®) x (9(10®)) matrices. These 
matrices are constructed by grouping together monomial terms that appear 
in eq. flTS]) and eq. flTS]) . RGIsearch performs this task by calculating these 
derivative monomial terms one-by-one and directly placing them into the 
matrix. To be able to do that, a std: :map is used to make the translation 
from a monomial term to a row in the matrix. 

Monomial terms are stored in a std: :vector of small integer numbers 
indicating the parameters and their powers. Because multiple numbers are 
required to fully specify a monomial term, and because containers such as 
vector require some memory overheac^l, the Cantor pairing function-. 

F : N X N —)■ N where P (ni, 77 - 2 ) = ^ (ui + 772 ) (771 + 772 + 1) + 772 , (17) 

is used to convert monomial terms into a single natural number, which is 
then mapped to a row in the matrix. To pair more than two numbers, the 
vector is extended with zeroes until its size is a power of two. The pairing 
function is then applied pairwise until a single number remains. The Cantor 
pairing function is injective, so in principle there is no danger of collision. 
However, the resulting number might overflow its designated 128 bit integer. 
Even if this happens, the probability of a collision is negligibly small|^ 


^Containers without overhead exist, but their basic functions such as sizeO, which is 
heavily used throughout the algorithm, have longer computation times, slowing down the 
algorithm significantly. 

^The calculation of this probability is comparable to the birthday problem. The 
probability of collision for a 128 bit integer and a matrix of typical size can be approximated 
to 1 - exp(-10iV2i29) 
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2.6. Markowitz Pivoting and Sparse Matrix Handling 

Next, RGIsearch solves the systems that were created in the above 
procedure. To hnd the nullspace of these systems, RGIsearch uses fraction 
free Gaussian elimination with Markowitz pivoting. An implicit pivoting step 
is performed before every elimination step to preserve sparsity. This pivoting 
is based on the heuristic Markowitz count j^: 

(g - 1) (g - 1 ), (18) 

where r* is the number of nonzero elements in row i and Cj the number 
of nonzero elements in column j. Before every elimination step, the nonzero 
element with the lowest Markowitz count is located and (implicitly) pivoted 
to the top-left. By performing this pivoting step, sparsity is preserved and 
and the number of required operations is reduced as much as possible. 

Note that the Markowitz count is zero for rows or columns with only a 
single nonzero element (singlets). These rows and columns are selected hrst, 
since they are very cheap to eliminate and produce no additional hll-in. The 
matrices that appear while computing polynomial invariants usually contain 
many row singlets, so the algorithm is specihcally tailored to use these. It is 
abstractly shown below.. The ith row of a matrix A is denoted with Aj, and 
the matrix is viewed as a collection of these rows since the order of the rows 
can freely be changed in nullspace computations. 

The reduction of a row simply means that all nonzero elements of the row 
are divided by their total greatest common divisor (gcd). If the algorithm 
terminates while A is not yet empty, then all variables have been eliminated 
and no nullspace exists. If A does empty, the nullspace can easily be solved 
by selecting a basis for the space of variables that have not been eliminated. 
The rest of the variables are then backwardly solved by the equations in B. 

For the method described in subsection 1.4, elimination of matrices of 
similar sizes containing univariate polynomials with integer coefficients has 
to be performed. The algorithm is very similar. The Markowitz count is 
extended to: 


(g - 1) (g - 1) {dij - 1), (19) 

where dij is the degree of the matrix element at row i and column 
j. Furthermore, integer coefficients are maintained through fraction-free 
Gaussian elimination using the subresultant polynomial remainder sequence 


Algorithm: Gaussian elimination with Markowitz pivoting 
Input : A sparse matrix A G 

Output: A fully eliminated matrix B G a vector of flags that 

indicate elimination of variables o G Z 2 
B i — Q i — 0, o — (0, •••,0) 

for all k G do 

if Ak is a row singlet then 
1^ Q Q U {k} 


while A is not empty or o still contains zeroes do 
i •(— 0, j ^ 0 
if Q 7 ^ 0 then 

i •(— Some element of Q 

Q Q — {i} 

if Ai is a row singlet then 

\_ j Column index of the nonzero element of Ai 

else 

Find the nonzero element Omn in A minimizing 

(r^ - 1) (Cn - 1) 
i m, j n 

if j 7^ 0 then 

Oj i — 1 

B ^ BU {Ai} 

for all I with aij 7 ^ 0 do 

Fraction-freely eliminate aij from Ai using row Ai 
Reduce Ap using the gcd of its elements 
if Aj is empty then 
lA^A- {A,} 

if Ai is a row singlet then 

|_ Q Q U {/} 

_ A ■<— A — {Aj} 
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algorithm outlined in . The eigenvalues are then computed by solving the 
polynomials that appear on the implicit diagonal of B. 

To improve speed and memory efficiency, RGIsearch uses the list of lists 
method to store matrices [sj. They are stored as a list of rows, containing 
all nonzero matrix elements which are tagged by their column index. While 
this method is not the most efficient of all sparse matrix storage techniques 
in terms of memory usage, it allows for very fast access to individual rows 
which is extremely important in the above algorithm. Additionally, insertion 
of new nonzero matrix elements is relatively cheap which is signihcant for 
the elimination steps. Fast access to columns is also important during the 
elimination steps. Therefore, a second matrix is stored as a list of columns 
instead of rows. This matrix is used purely as an indexing structure and thus 
only stores the row indices of nonzero elements. Both structures are updated 
simultaneously during the elimination procedure. 

2.1. Two-loop invariants 

RGIsearch allows for computation of regular polynomial invariants for 
two-loop /9-functions. The general form of a two-loop /3-function is: 

/3 (xi) = (xi) + • (20) 

Therefore, a two-loop invariant looks like: 

Taking the derivative and letting denote the part of the derivative 
involving the condition for invariance for the one-loop and two-loop 

terms separate: 


= 0 and + 4^^ = 0. (22) 

The term 4^^ is formally of three-loop order. These conditions are 
converted into a linear system in a very similar fashion to the one-loop case. 

3. Usage 

RGIsearch is available from a public repository [^. The latest version of 
RGIsearch can be acquired with: 

$ git clone https://github.com/rbvh/RGIsearch.git 
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3.1. Compilation 
To compile, do: 

$ make 

This will generate all binary files and link them into an executable. This 
executable can then be run with; 

$ ./RGIsearch 

If any /9-functions, settings or any of the sources are changed, make 
should update the binaries. If for some reason it does not, or it is otherwise 
required, one can do: 

$ make clean 
$ make 


3.2. Settings 

The program settings can be found in the hie settings . h . 

const int DIM_SEARCH_PARAMETER 

This constant controls the range of dimensionalities the algorithm 
searches. Depending on the dimensionalities that are found and the value 
of this constant, RGIsearch calculates a suitable selection of dimensionality 
vectors to search through. 

const int MAX_TERM 

This constant controls the number of different parameters allowed in 
a monomial term during the execution of the algorithms for polynomial 
invariants. For instance, if it is set to 2, terms like xy"^ can occur, but xyz is 
forbidden. Ghanging it can have huge inhuence of the required computation 
times and memory. 

const int FILTER_THRESHOLD 

This constant controls the hltering mechanism. If it is set to a higher 
value, the hltering algorithm will perform a more elaborate attempt to hud 
all possible products of previously found invariants to compare against any 
newly found invariants. The default value should almost always be sufficient. 

bool INCLUDE_TW0_L00P 

Set to true to search for invariants at two-loop level. 
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bool FACTORIZE 


Set to true to include the algorithm that searches for factorized polynomial 
invariants. The code does currently not support factorization for two-loop 
level searches. 

bool REPORT 

Set to true to make the program report its activities in more detail. 

3.3. (3-Functions 

The /3-functions can be found in equations.cpp . 

To dehne the /9-functions of a theory, the parameters of the theory hrst 
need to be dehned. This is done with: 

Param newPar("parName" , parSize); 

parSize is the size of the parameter (1 for a scalar, n for anxn matrix). 
If a parameter is complex, its daggered counterpart must be defined. This is 
done with: 

Param newParDagger = dagger(newPar); 

The Param class has several methods that simplify matrix parameters. 
They include: 

1. botRightO sets all elements except for the bottom right component 
equal to zero. 

2 . diagO sets all offdiagonal elements equal to zero. 

3. botRightDiagO sets all offdiagonal terms equal to zero, and makes all 
components except for the bottom right degenerate. 

After dehning the parameters, one can dehne the /9-functions. These are 
stored as a vector of the objects BetaFunc: 

vector<BetaFunc> nameOfBetaFuncs 

The /3-function of a parameter can then be dehned as: 

BetaFunc bNewPar(newPar); 
bNewPar = <polynomial>; 
nameOfBetaFuncs.push_back(bNewPar); 

The polynomial can be constructed using regular arithmetic involving 
the parameters of the theory. Irrational numbers a/b can be represented as 
ir(a,b). A trace function is available as Tr(). For the complex conjugate 
parameters, their /3-functions must be included as: 
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BetaFunc bNewParDagger = Conjugate(bNewPar); 
nameOfBetaFuncs.push_back(bNewParDagger); 

RGIsearch needs separate vectors of /3-functions for one-loop and two- 
loop, which must be in the same order. Finally, the algorithm can be initiated 
by calling: 

findinvariants(lloopBetaFuncs, 21oopBetaFuncs); 

3.4- Test Program 

As an example, we consider a simple toy system of /3-functions for three 
variables x, y, z: 



(23a) 


(23b) 


(23c) 


This system has two invariants at two-loop level. To hnd them, the 
appropriate equations.cpp is shown below. 

Using the default settings and enabling two-loop calculations, the output 

is: 


—X"2 -F 2y''2 — z -F (l/16pi''2){x -F y} 
—2y''3 — 2xz -F (1 /16 pi" 2) {x"2 -F 3z} 


These can easily be verihed to be invariants of the above system. 
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Test case equations, cpp file 

#include "src/common.h" 

#include "src/invariants.h" 

int mainO 

{ 

Param x("x", 1); 

Param y("y", 1); 

Param z("z", 1); 

vector<BetaFunc> test_l; 
vector<BetaFunc> test_2; 

BetaFunc bx_l(x); 

BetaFunc bx_2(x); 
bx_l = -4*x*y - 3*y*y; 
bx_2 = -l*x + 6*y; 
test_l.push_back(bx_l); 
test_2.push_back(bx_2); 

BetaFunc by_l(y); 

BetaFunc by_2(y); 
by_l = -2*x*x + l*z; 
by_2 = 3*x; 

test_l.push_back(by_l); 
test_2.push_back(by_2); 

BetaFunc bz_l(z); 

BetaFunc bz_2(z); 
bz_l = 6*x*y*y + 4*y*z; 
bz_2 = -4*x*y - 3*y*y + l*z; 
test_l.push_back(bz_l); 
test_2.push_back(bz_2); 

findinvariants(test_l, test_2); 

} 
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